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ABSTRACT 

We present 2D axisymmetric hydrodynamic simulations investigating the long term 
effect of FR II radio galaxies on the metal distribution of the surrounding intra-cluster 
medium (ICM). A light jet is injected into a cooling flow atmosphere for 10 — 30 Myr. 
We then follow the subsequent evolution for 3 Gyr on a spherical grid spanning 3 Mpc 
in radius. A series of passive tracer particles were placed in an annulus about the cluster 
core to simulate metal carrying clouds in order to calculate the metallicity (Z) as a 
function of time and radial distance from the cluster centre. The jet has a significant 
effect on the ICM over the entire 3 Gyr period. By the end of the simulations, the jets 
produced metallicities of ~ 10% of the initial metallicity of the cluster core throughout 
much of the cluster. The jets transport the metals not only in mixing regions, but also 
through upwelling ICM behind the jet, enriching the cluster over both long and short 
distances. 
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1 INTRODUCTION 

The intra-cluster gas has a metalli city of a thir d the so- 
lar v alue (Zq) on average (e.g. iMolendil l20ui iRenzinil 
l2004h . with values of « 0.2 Zq up to the megaparsec scale. 



There is also evidence that the metals were expelled from 
th e most massive galaxies in a cl uster at redshifts z > 
1 ([Mushotzkv &; Loew enstein 1997) rather than being re- 



moved by ram pressure (Renzini 2004). The energy of the 



supernovae that produced these metals is hardly enough to 
expel the metals from the galaxies. An obvious alternative 
are radio sources produced by jets from the active galac- 
tic nuclei of giant ellipticals. Weak jets, as often found in 
nearby clusters , have been found in efficient in this respect 
( Briiggenl 120021: lOmma et al.1 l2004h . FR II radio sources 
Fanaroff fe Rilevlll974h are certainly powerful enough to do 



the job, especially at high redshift (e.g. ICarilli et al.l 1200 ll ). 
However, the actual mechanism by which the jets would do 
this is so far unclear. 

It is becoming increasingly clear, from both obser- 
vational and theoretical work that radio sources have a 
significant impact on their environment. The y push aside 



120031 : iKrausd 2003). After the central engine has shut 
down, the low entropy plasma that fills these cavi- 
ties continues to rise through the surrounding atmo- 
spheres due to buo y ancy and excess momentum (e.g . 
I Churazov et al.1 l200ll : iRevnolds et~aT1 12002b iKaiserl 120031 : 
I Kaiser et al. 20051 ). This may lead to megaparsec scale con- 
vection, as demonstrated recently in 3D hydrodynamic sim- 
ulations (jBasson fe Alexander] [20031 ). 

In the following we investigate if this convection is able 
to eject the metals fr om a giant elliptical galaxy. We follow 
a similar approach to iBasson &; Alexander! (|2003l ) , with the 
primary difference being the inclusion of tracer particles, to 
simulate metal containing clouds. 

In contrast to previous work (|Briiggenl 120021 : 
lOmma et al.l |2004) our approach is tailored for high 
power FR II sources. 
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The associate d heating is relevant to the cooling flow 
problem (e.g. iKaiser &; Alexander! Il999l : iKaiser &; Binnevl 
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2 THE SIMULATIONS 

The hydrodynamic simulations were constructed to model 
an FR II radio jet propagating through a cooling clus- 
ter atmosphere. The conservation equations for mass, mo- 
mentum and energy were solved axisymmetrically, in two- 
dimensional spheri cal polar coordina t es (r, 0), employin g 
the FLASH code (|Frvxell et al.1 l2000l : ICalder et al.fl2002h . 
FLASH is a conservative, adaptive-mesh hydrodynamics 
code using third order interpolations. The simulated domain 
is a sphere with radius 3 Mpc, with a 20kpc (= Tinner) ra- 
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Table 1. Simulation parameters 



1 x10^ 4 2x10^' 



3x1 24 4x1 2 ' 
x (cm) 



5x10 24 6x10 2 ' 



time - 3000.047 Myrs 

number of blocks = 2189, AMR levels = 7 

Figure 1. Block distribution for a typical run. Each block is 
resolved by 8 x 8 cells. The background displays the density, where 
the highest values are in the centre (2 x 10 -24 kgm -3 , blue) and 
the lowest ones in the remains of the rising buoyant cocoon (2 X 
10 -28 kgm -3 , red). Only the central 2x4Mpc of the 2D-spherical 
grid are shown. 



dius spherical cutout at the centre, with outflow boundary 
conditions (zero gradient, velocity unrestricted) set for when 
the jet is not active, allowing material to pass across it. The 
outer boundary at radius 3 Mpc is set to outflow also, so 
material can flow in from outside the simulation, whilst the 
= and = tt boundaries are set to reflective to account 
for the rotational symmetry. The grid uses an adaptive mesh 
with a maximum of seven refinement levels to give effective 
maximum resolution of 6kpc and 0.1 kpc in the radial and 
angular directions, respectively; this gives up to 512 cells in 
both directions. The code employs square blocks of eight by 
eight cells. A typical distribution is shown in Fig. [1] 

The intra-cluster medium (ICM) was modelled as a 
monatomic, isothermal (10 8 K) King profile atmosphere, 
with a density profile given by 

-3/3/2 
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Run A 
Run B 


10 

30 


3000 
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3 x 10- 26 
1 x 10- 26 
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4.7 x 10 39 
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a Jet's active time. 

^ Total simulation time. 

c Mass density in the beam at inlet. 

^ Total energy released by both jets. 

e Luminosity per jet, L jet = 7rr f nn er^jet' y fet( 1 - cos 



where /3 was chosen to be 2/3 and the core radius (oo 
and core density (po) were 100 kpc and 10" 23 kg m 3 re 
spectively. A gravitational acceleration of 



9r 



3f3k B T r 
/im H a 2 [l + (r/a ) 2 ] 



(2) 



towards the centre is used to keep the uncooled ICM in 
hydrostatic equilibrium, where \i — 1/2, /3 = 2/3 and ttih is 
the mass of h ydrogen. A Bremsstrahlung cooling term (e.g. 
lLongaiil li994) was included of the form 



1.722 x io~ 4 Wtwi 



(3) 



P\cm = po (l + (r/a ) 



(1) 



dE 
dt 

assuming the ICM is comprised of ionised hydrogen and 
where n is the particle number density. 

This setup was run for 200 Myr to allow an initial cool- 
ing flow to establish. To follow the distribution of metals 
~ 10 4 passive tracer particles (chosen to obtain good statis- 
tics) were positioned between radii of 40 to 150 kpc from the 
cluster core so as to represent constant initial metallicity in 
the shell. The particles are introduced after 200 Myr, once 
the cooling flow is established - subsequently the particles 
move with the local velocity of the ICM. 

Once the initial cooling flow has been established and 
particles added, the radio source is introduced by injecting 
low density material across the inner (r = 20 kpc) boundary 
as two anti-parallel jets in pressure equilibrium. Each jet has 
a half-opening angle of # op en = 0.5(= 29°), and a velocity of 
c/3. The jets are therefore supersonic with respect to both 
the external and internal sound speed. 

We present two simulations. They have the same to- 
tal energy but differ in the instantaneous kinetic jet power, 
modelled by a trade-off between jet density and active time. 
Important parameters are shown in Table [T] for each run. 



3 RESULTS 

Logarithmic density maps, including particle distribution, 
are shown in Figs. [2] and |3] The general dyn a mics is very 
similar to that found by iBasson Sz Alexander! (|2003h . Dur- 
ing the active phase, the jet develops the shock and cocoon 
structure that is typical for such simulations. After the jet 
activity has ceased, the centre quickly refills, and a connec- 
tive flow establishes. The ICM flows outward near the jet 
axis, and comes back near the equatorial regions. 
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Run A: 1000 Myr 
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Run A: 2000 Myr 
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Run A: 3000 Myr 
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Figure 2. Logarithmic density maps for run A at four times. Positions of the passive particles are overlaid on the right hand side for 
each time. 



3.1 Particle transport 

The particles are not picked up by the bow shock, as one 
might naively have expected. This is easily explained by 
realising that the front speed is generally larger than the 
gas velocities. The front therefore quickly passes, whereas 
gas and particles do not change their position significantly. 
The particles generally honour the contact surface, and are 
therefore displaced by the jet's cocoon. This can be observed 
early in the evolution and also in later phases (Figs. [2] and [3} 



- this demonstrates the accuracy of the particle transport. 
The particles do enter intermediate density regions, where 
the jet plasma has mixed with the surrounding gas. As 
the simulation proceeds, instabilities cause the cocoon to 
loose its integrity and to "dissolve" into a "bubbly" struc- 
ture. The tracer particles are found in the ICM gas be- 
tween the bubbles, and as the turbulent structure rises 
buoyantly the tracer particles are dragged from the clus- 
ter core. At this time we also see a large-scale convective 
now established which is very similar to the one found by 
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Figure 3. Same as Fig. [2] but for run B. 



iBasson fe Alexander] (|2003h . The particles are effectively en- 
trained in the flow. In both simulations, the particles easily 
reach megaparsec scales at late times. 

The particles show an artificial feature in the plane of 
symmetry. In the simulations presented here, there is very 
little asymmetry in the initial conditions. We also performed 
a simulation in which we introduced a significant asymmetry 
between the two jets, but which is in other respects identical 
to run A - the feature in the symmetry plane is then strongly 
suppressed. About 10% of the particles are in this feature, 



and it is significant in the region out to « 100 kpc, only. 
About 1% of the particles are lost across the inner radial 
boundary during the simulation. 

Figure U shows the final versus initial radii for all the 
surviving particles for runs A and B. In both cases, most 
of the particles have been displaced by about a megaparsec. 
Their final position is independent of their initial position. 
We conclude that our results are likely to be insensitive to 
the precise initial distribution chosen for the tracer particles. 
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Figure 4. Final versus initial particle radii for run A (left) and run B (right). 



3.2 Metallicities 

The tracer particles can be thought to represent a certain 
mass of metals, which can be assigned to each particle indi- 
vidually. In the following we carry out a metallicity analysis 
with the initial 3D metallicity distribution normalised to be 
a constant in space. Since the initial particle distribution is 
homogeneous, the metal mass for each particle of the initial 
distribution is given by: 

^particle OC V (4) 

1 + (r/ao) 

This mass is a property of the particle and is retained as it 
moves during the simulation. In order to derive the metal 
mass in a given volume, we sum the contributions from every 
particle. For the metallicity in a 3D radial shell, this gives: 

Z ^ Sshell volume particle , . 

fv, 11 i P dV 

J shell volume ' 

Figure [5] shows the time evolution of the metallicities 
relative to the initial core metallicity against radius for both 
runs. Due to the reduction of the core gas density, the metal- 
licity first increases at small radii by a factor of a few, then 
decreases as the particles move outwards. In run A, the cen- 
tral metallicity drops below its initial value, whereas in run B 
it stays above its initial value throughout the simulation. 
In both simulations, the metallicity reaches « 5-10% of the 
initial core metallicity on the megaparsec scale. The high 
power jet (run A) distributes the metals over a somewhat 
larger range out to larger radii than the lower power jet. The 
final metallicity distribution drops by a factor of ten out to 
one megaparsec for Run A, and is almost constant between 
0.4 and 0.8 Mpc for Run B. 



4 DISCUSSION 

We have presented simulations of powerful jets which com- 
prise an active phase as well as a phase long after the activity 
of the central source has ceased. Metal carrying clouds are 
simulated by a set of passive tracer particles. For the pur- 
pose of the present discussion, it is not relevant if the clouds 



remain intact or are instead dispersed, since we are only in- 
terested in the metallicity averaged over significant volumes. 
However, the cloud mass may be important. The energy to 
lift the clouds out of the cluster core comes from the kinetic 
energy of the convective flow. Hence the mass of the clouds 
must be small compared to the total gas mass for our use 
of passive tracer particles to accurately represent the evo- 
lution of metallicity - this should be the case for realistic 
situations. 

The particles are dragged out to the megaparsec scale 
by the large-scale convective flow and between the radio 
bubbles. Towards the end of the simulations they occupy 
a linearly shaped region. The rising timescale (Gyr) is of 
the order of the crossing timescales for galaxies in the clus- 
ter. These will disperse the metals further. The same would 
be true for subcluster mergers. The exact metallicity distri- 
bution resulting from our simulations might therefore not 
be directly observable. 

Quantitatively, the two runs differ slightly. Both simu- 
lated jets release the same total amount of energy. Run A, 
having a higher power but shorter active time, expels 
nearly all the tracer particles from the host galaxy, and 
distributes them to somewhat larger radii, with generally 
slightly smaller metallicity than in Run B. 

This suggests, that high power sources might be more 
effective in the metal enrichment of clusters. Such sources 
are predominantly found at high redshift consistent with 
the observation that metal enrichment should have taken 
place at z > 1 ([Mushotzkv &; Loewensteinlll997h . The cen- 
tral metallicity gradient is quite different in run A compared 
to run B: in run A the central metallicity profile is quite flat, 
while in run B a steep central metallicity profile persists. In- 
terestingly, observed cooling flow clusters show steep central 
metallicity gradients, while in non-cooling flow clusters it is 
rather flat. This has b een explained by recent clust er merg- 
ers in the latter cases (|De Grandi fe Molendi 2001) and the 
absence thereof in the former. Our results suggest that only 
a fairly high power radio source flattens the metallicity gra- 
dients. Note that the radio sources at the centre of cooling 
flows in the local universe are almost exclusively of FR I 
type, i.e. they would have much lower power than the jets 
simulated here. Hence, we would not expect that they would 
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Figure 5. Derived metallicities for different simulation times for run A (left) and run B (right). The metallicities are normalised to the 
initial core metallicity. The symmetry feature (see text) shows up in a narrow range around 100 kpc, for late times. 



erase central metallicity grad ient s. This is con s istent with 
the findings of lBriiggenl (|2002h and lOmma et all (|2004h , who 
simulate the metal evolution for low-power, subsonic flows. 
If a cluster would have produced a radio source with a power 
of more than 5 x 10 39 W with a lifetime of about 10 Myr, 
within the preceding 2-3 Gyr, which is our simulation time, 
we would expect an effect on the central metallicity gradi- 
ent. The timescale is comparable to the time needed by the 
host galax ies to build up the metallicity gradie nts via stellar 
evolution ([David et al.ll200ll : IWise et al.l l2004). so the exact 
gradient would depend on the details. 

Overall, the metallicities come close to 10% of the initial 
core metallicity. Therefore, depending on the metallicities 
in the host galaxy, one outburst might already produce a 
significant fraction of the metals observed in present day 
clusters. 

Although the presented simulations are 2D axisymmet- 
ric, we expect the simulations to recover the basic trends 
from full 3D simulations, since computations similar to 
those presented here have already been performed in 3D 
(Bas son fe Alexand er 2003), and a very similar overall flow 
structure has been reported. A significant magnetic field 
would change the details of the mixing between the jet 
plasma and the cluster gas but could hardly influence the 
large scale flow structure. 



5 CONCLUSIONS 

We have shown that the remnants of powerful FR II jets 
can efficiently drag passive tracer particles out of the host 
galaxy. These particles are a good representation of metal 
containing clouds as long as the total mass of the clouds is 
small compared to the gas mass of the galaxy cluster. The 
resulting metallicity distribution was very shallow with gas 
on the megaparsec scale reaching metallicities of up to 10% 
the initial metallicity of the cluster core. The angular metal- 
licity distribution is peaked towards the former jet axis. Ne- 
glected processes in the cluster are expected to produce a 
more isotropic distribution. Observations generally show a 
homogeneous metallicity distribution, with values of about 
Z©/3, constant out to large radii. The metals are believed to 



originate predominantly from the largest cluster members, 
having been expelled at redshift z > 1. Powerful jets are 
detected out to much higher redshift s, and are believed to 
originate in the most massive galaxies at their redshift. As 
we have shown here, their remnants are able to expel met- 
als from their host galaxy and distribute them throughout 
the galaxy cluster. This provides additional evidence that 
FR II radio jets may contribute significantly to the metal 
transport in galaxy clusters. 
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